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Long baseline laser interferometers used for gravitational wave detection have 
proven to be very complicated to control. In order to have sufficient sensitivity 
to astrophysical gravitational waves, a set of multiple coupled optical cavities 
comprising the interferometer must be brought into resonance with the laser 
field. A set of multi-input, multi-output servos then lock these cavities into 
place via feedback control. This procedure, known as lock acquisition, has 
proven to be a vexing problem and has reduced greatly the reliability and duty 
factor of the past generation of laser interferometers. In this article, we de- 
scribe a technique for bringing the interferometer from an uncontrolled state 
into resonance by using harmonically related external fields to provide a deter- 
ministic hierarchical control. This technique reduces the effect of the external 
seismic disturbances by four orders of magnitude and promises to greatly en- 
hance the stability and reliability of the current generation of gravitational 
wave detector. The possibility for using multi-color techniques to overcome 
current quantum and thermal noise limits is also discussed. 
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1. Introduction 

Gravitational waves promise to reveal new information 
about the bulk motions of massive compact objects in 
the universe. In this decade kilometer-scale interferom- 
eters, such as LIGO [1, 2], Virgo [3], GEO600 [4], and 
KAGRA [5], are expected to make the first direct de- 
tection of gravitational waves in the 10- 10k Hz band. 
The Advanced LIGO (aLIGO) project [6] is a signifi- 
cant upgrade of the initial LIGO interferometers, includ- 
ing more sophisticated vibration isolation, a factor of ten 
higher laser power, larger test masses, and a more versa- 
tile optical readout, among other improvements. These 
improvements should lead to a factor of 10 sensitivity im- 
provement across the entire detection band, resulting in 
a factor of 1000 in increased probed volume of space. 

The problem of moving an interferometer from its ini- 
tial uncontrolled state (where the suspended mirrors are 
swinging freely) to the final operating state (where all 
cavity lengths are interferometrically controlled) is re- 
ferred to as "lock acquisition" . For a single Fabry-Perot 
cavity or a simple Michelson interferometer, the prob- 
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lem is relatively straightforward: typically the locking 
servo is engaged and the mirrors are moved until the 
servo forces the interferometer into the desired operating 
point. With a more complicated configuration utilizing 
multiple coupled cavities (e.g. aLIGO) there are no re- 
liable optical signals providing cavity length information 
until all cavities are simultaneously resonant. This prob- 
lem is compounded by seismic-induced residual mirror 
motions which, even with the advanced seismic isolation 
systems used in aLIGO, are expected to be ~ 10~^ m 
below 1 Hz [7] . Waiting for full acquisition to happen by 
chance is an exercise in futility. 

Over the years various techniques have been developed 
to address the lock acquisition problem. Algorithms have 
been developed that use digital controls, clever sequenc- 
ing of feedback loops and mixing of interferometric sig- 
nals to reduce the waiting time [8-10]. Unfortunately, 
these are still insufficient for the needs of aLIGO. In- 
stead, aLIGO will be using an auxiliary arm length sta- 
bilization system to robustly bring the long arm cavities 
to a stable operating point at the cavity resonance, inde- 
pendent of the rest of the interferometer. This technique 
uses frequency-doubled auxiliary lasers, phase locked to 
the main interferometer laser source, to serve as an in- 
dependent sensor for the cavity lengths. The primary 
motivation of the experiment described in this article is 
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to demonstrate that the arm cavity length can be in- 
dependently controlled by auxiliary locking to within a 
small fraction of the linewidth of the arm cavity at the 
primary laser frequency. 

Beyond the practical application of aLIGO arm length 
stabilization, this type of "multi-color metrology" can 
allow us to make better measurements of cavity proper- 
ties and noise fluctuations. This may one day make it 
possible to sense inherent noise sources and reduce their 
metrological effects. 

Arm length stabilization with frequency-doubled auxil- 
iary lasers has been demonstrated by Mullavey et. al [11] 
in shorter-baseline cavities. Other methods for indepen- 
dent arm cavity stabilization have included digital inter- 
ferometry [12, 13], and suspension point (or suspension 
platform) interferometry [14, 15]. 

This article describes a prototype experiment of cavity 
length stabilization using multiple laser wavelengths, and 
its implications for future interferometers. The experi- 
ment has been performed on a 40-meter-long suspended 
Fabry-Perot cavity. Possible noise sources and their con- 
tributions are discussed. 

2. Experimental Setup 

This experiment was conducted on the Caltech-LIGO 
40m prototype interferometer [16-18]. This prototype is 
used to develop interferometer topologies for future grav- 
itational wave detectors. Currently the configuration is 
similar to that of aLIGO (dual-recycled Michelson with 
40m-long Fabry-Perot arm cavities), and is being used to 
prototype aLIGO readout and control schemes. "Dual- 
recycling" refers to the use of both power and signal recy- 
cling mirrors at the input and output ports respectively of 
the Michelson interferometer, recycling light power that 
would otherwise escape through those ports [19, 20]. All 
of the main interferometer optics are suspended as single- 
stage pendula with a length of 25 cm {fp = 1 Hz). The 
entire interferometer is housed in an ultra-high vacuum 
envelope. For this experiment just a single, suspended, 
40m-long Fabry-Perot arm cavity of the full interferom- 
eter is used. The rest of the interferometer optics are 
misaligned so as to not affect the measurements. 

Figure 1 shows a schematic diagram of the experimen- 
tal setup. A 1064 nm beam (red in diagram) from the 
main interferometer pre-stabilized laser (PSL) is injected 
into the arm cavity through the input mirror (ITM) from 
the interferometer vertex. A second, 532 nm beam (green 
in diagram) from an auxiliary (AUX) laser is injected 
through the cavity end mirror (ETM). Whereas the PSL 
beam circulates through the full interferometer under 
normal operating conditions, the AUX laser beam res- 
onates only in the single arm cavity and is extracted be- 
fore interacting with the rest of the interferometer. 

The cavity mirrors are dichroic and highly reflective 
at both wavelengths. For the 1064 nm beam, the cav- 
ity is overcoupled and reflects most of the light back to- 
wards the interferometer vertex. For the 532 nm beam 
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Table 1. Cavity properties of the arm cavity, at 1064 nm 
as seen from the vertex, and 532 nm as seen from cavity 
end. 



the cavity has a much higher transmissivity and some of 
the light is transmitted through the cavity at the ITM, 
through the interferometer vertex and extracted from the 
vacuum. The mirror and cavity properties for both wave- 
lengths are shown in Table 1. 

A. PSL light source 

The 1064 nm light source is a 2 W Innolight NPRO. The 
light is spatially filtered by a ~ 20 cm ring cavity with a 
2 MHz bandwidth that also provides passive filtering 
of the laser noise at RF frequencies. The PSL frequency 
is locked to an in-vacuum, suspended, critically-coupled, 
triangular mode cleaner (MC) cavity, which conditions 
the beam by suppressing excess frequency noise and re- 
jecting higher-order spatial modes. The laser is locked 
to the MC via the standard Pound-Drever-Hall (PDH) 
method [21] with a ~ 130 kHz bandwidth. The power is 
adjusted to allow approximately 25 mW of 1064 nm laser 
light to be incident on the cavity under test. 

B. A UX light source and frequency doubling 

The AUX beam comes from a frequency-doubled 700 mW 
JDSU NPRO-126N. The frequency doubling is achieved 
via second harmonic generation (SHG) in a PPKTP crys- 
tal [22]. The conversion efficiency is ~1%/W, and with 
other input losses we end up with 1.2 mW of 532 nm 
light incident on the ETM. 

C. Dichroic mirror coatings 

The mirrors forming the Fabry-Perot arm cavities have 
custom coatings to provide reflectivity at both 1064 and 
532 nm. Figure 2 shows the calculated coating reflectiv- 
ity as a function of wavelength for the ITM (the ETM 
shows a similar profile). The layer structure is a partic- 
ular aperiodic design chosen to minimize the influence of 
various types of thermal noise on the reflected phase of 
the laser field [23]. 
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Fig. 1. Experimental setup. Red lines indicate the path of the 1064 nm PSL beam, and green lines indicate the path 
of the 532 nm AUX beam. The colored regions correspond to logical sections of the control and readout, described in 
more detail in Section 3. Optics in shadow are part of the larger interferometer not used in this experiment. 




Wavelength [nm] 

Fig. 2. Spectral reflectivity of one of the dichroic cavity 
mirrors. 



D. Sensing, acquisition and control 

Initially, the AUX laser frequency is locked to the arm 
cavity length via a standard PDH locking scheme. The 
AUX laser is locked to the cavity, rather than vice versa, 
because the laser frequency actuator has much greater 
bandwidth than the cavity mirror displacement actua- 
tors. Phase modulation sidebands at 217 kHz are in- 
troduced on the AUX beam by directly driving the laser 
frequency actuator. This modulation frequency is chosen 
to minimize the ratio of amplitude to phase modulation. 



The green light reflected from the ETM is used for the 
PDH lock. 

Once the AUX laser is locked to the cavity, the AUX 
beam transmits through the ITM and is extracted from 
the vacuum system. The extracted AUX beam is used 
in a heterodyne measurement with a frequency-doubled 
sample of the PSL beam. The frequency of the beat 
note between the AUX laser and the frequency-doubled 
PSL is measured by the a delay-line frequency discrim- 
inator (DFD) (see section 2E). The DFD has "coarse" 
and "fine" paths which provide different dynamic ranges. 
These outputs are the primary error signal for the cav- 
ity control. They are digitized and a control signal is 
generated with a digital feedback control system. 

Figure 3 shows the control sequence and hand-off be- 
tween the coarse and fine discriminator paths. Since seis- 
mic noise acting on the length of the MC and arm cavi- 
ties causes the beat note to fluctuate by about 10 MHz, 
the large range coarse path is used to engage feedback 
smoothly. After length control is achieved an artificial 
offset is introduced in the discriminator signal to sweep 
the length of the arm cavity until the length meets the 
resonance condition for the PSL beam. This ability to 
tune the cavity length directly is the key to the use of this 
technique as a lock acquisition tool for Advanced LIGO. 

In the end, control is passed to the fine discrimina- 
tor by digitally fading over between the coarse and fine 
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Fig. 3. Sequence of the arm length control as a function of time. The intra-cavity power and detuning are for the 
1064 nm PSL beam. 



signals. The fine discriminator has a higher signal-to- 
noise by a factor of the extra delay. At this stage the 
arm length can be tuned more precisely so that the main 
laser fully resonates in the arm. The final steady-state 
control used during the measurement is described in de- 
tail in section 3. A PDH error signal derived from the 
PSL beam refiected off of the ITM provides an out-of- 
loop measure of the residual cavity displacement. 

E. Delay-line frequency discriminator 

The delay-line frequency discriminator (DFD), which is 
used to measure the frequency of the beat note between 
the PSL and AUX beams, works by mixing an input RF 
signal with a delayed version of itself. For a given delay in 
the delay line, the mixer output voltage is a periodic 
function of the input frequency, /. In the small frequency 
limit (/ <C l/D) the output is directly proportional to 
the input frequency: 



mixer low pass preamp 



V oc Df. 



(1) 



Figure 4 shows a schematic of the DFD circuit. The 
signal from a broadband RF photodetector first passes 
through a comparator that turns the signal into a square 
wave. This helps reduce noise associated with small 
amplitude fiuctuations of the input signal. This signal 
is amplified and split into two discriminator paths: a 
coarse path with a delay of 7.3 ns and frequency range 
of 34 MHz, and a fine path with a delay of 270 ns and a 
range of 3.6 MHz. The mixer outputs, with signals given 
by equation 1, are filtered and digitized and used as er- 
ror signals for the cavity length control servo. The coarse 
path, which has a larger bandwidth, is used during lock 
acquisition, whereas the lower-noise fine path is used to 
achieve best performance in the steady state. 
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Fig. 4. Delay-line frequency discriminator. Compo- 
nents used: comparator: Analog Devices AD9696; RF 
amplifier: Mini-Circuits ZHL-IA; power-splitter: Mini- 
Circuits ZBSC-413+; delay line cables: RG-58 C/U; 
mixer: Mini-Circuits ZP-3+; low-pass filters: Mini- 
Circuits BLP-1.9; pre-amplifier: Stanford Research 
SR560. 



In order to confine the cavity length to the within the 
line width of the PSL a residual fiuctuation level of 10 pm 
RMS must be achieved. This means that the frequency 
noise of the fine path needs to be less than 7.4 Hz/V^Hz 
in the control bandwidth. The noise of the two paths 
are currently limited by the active readout electronics at 
an estimated level of 2.0 Hz/\/Hz and 0.1 Hz/V^Iz for 
the coarse and fine paths respectively. They therefore 
reasonably meet the required frequency stability. 

A delay-line design was used, rather than a phase- 
locked loop (PLL) design [24] because the frequency 
range of DFDs are relatively easy to tune and can be 
adjusted to give a large frequency range (e.g. the coarse 
channel). DFDs also don't require any active feedback 
loops, which complicate PLLs. Alternatively, a combi- 
nation of a large range DFD and a smaller range PLL 
could be a possible solution depending on the required 
frequency range and noise. 
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3. Control Model 

In this section we present a model of the control system 
used in the experiment. The control system, shown in 
Figure 5, can be broken into five parts, each described in 
the sections below. The model includes injection points 
for various noise sources that might affect overall perfor- 
mance, discussed in detail in section 4. 

A. AUX-cavity loop 

The first logical control loop is the PDH lock of the 
frequency-doubled AUX laser to cavity (green block in 
Figures 1 and 5). This loop suppresses the frequency 
noise of the AUX laser and allows its frequency to follow 
the motion of the cavity length. The control bandwidth 
of this loop is 30 kHz, limited by the laser cavity PZT 
frequency actuator. 

Since the AUX laser is locked to the cavity, information 
about the length fluctuation of the cavity is encoded in 
the frequency of the AUX laser light transmitted through 
the cavity ITM. The relationship between cavity length 
and laser frequency is given by: 



dL 



du 



(2) 



where L is the cavity length and u is the frequency of 
the laser resonating in the cavity. The cavity frequency 
response is a function of the cavity finesse and can be 
approximated as a single-pole low pass filter: 



1 



Kf/fcY 



(3) 



where / is the frequency of the signal and fc is the cav- 
ity pole frequency, which for the 532 nm AUX beam is 
18 kHz. 

The photo detection and mixing process that produces 
the PDH error signal has an overall flat V/Hz conversion 
factor given by I^aux- The servo filter, which is tuned 
to provide stable and robust locking, has a frequency 
response of Faux- Finally, the laser PZT frequency ac- 
tuator has a response of Aaux- 

The noise sources associated with this loop are the 
AUX laser frequency noise, second-harmonic generation 
noise, shot noise at the detector, and electronics noise 
of the readout electronics. These noise sources will be 
discussed in detail in section 4. 

B. PSL-MC loop 

The PSL-MC loop (yellow in Figures 1 and 5) describes 
the lock of the PSL frequency to the mode cleaner length. 
The control bandwidth of this loop is 130 kHz. As with 
the AUX laser cavity loop, the MC cavity has a single- 
pole frequency response given by Cpsl- ^psl is the re- 
sponse of the PDH sensing, Fpsl is the response of the 
servo filter, and ApsL is the response of the PSL com- 
pound frequency actuator, which includes the laser crys- 
tal temperature actuator, laser cavity PZT, and an EOM 
phase modulator. 



The PSL light transmitted through the MC and inci- 
dent on the arm cavity under test is the reference for the 
performance of the cavity stabilization system. However, 
finite gain in the PSL-MC loop can potentially lead to 
PSL frequency noise coupling into the arm stabilization 
loop. Despite this, the model shows that the suppression 
ratio from PSL frequency noise to residual displacement 
is more than 10^ at 100 Hz, so we can safely neglect the 
PSL frequency noise. 

C. Beat note detection 

The beat note detection block (blue in Figures 1 and 5) 
measures the frequency difference between the AUX and 
frequency-doubled PSL beams. 

The transmitted AUX light is a combination of the 
suppressed AUX laser frequency noise and any external 
displacement noises in the cavity that modulate the op- 
tical phase of the laser resonating in the cavity. These 
external disturbances result in phase noise that is con- 
verted to frequency noise by multiplying by the complex 
frequency, if. This signal is then low-pass filtered by the 
cavity pole, Caux, with a resultant transfer function of: 



AUX 
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i(///c) 
l + i(///c 



(4) 



where again = 18 kHz. I^dfd represents the flat 
Hz^V conversion of the full beat detection process, in- 
cluding the gain of the RF photo detector and the delay- 
line frequency discriminator. 

Noise in the frequency doubling process of the PSL 
should be at a similar level to that in the AUX-cavity 
loop. Laser shot and photo-detection noises are also sim- 
ilar to those in AUX-cavity loop. The electronics noise 
in this case is from the frequency discriminator (DFD). 
Finally, there is also additional noise from the analog-to- 
digital conversion (ADC) process. 

D. Arm length servo 

The output of the beat note detection process is the er- 
ror signal for the arm cavity length control servo loop 
(pink in Figures 1 and 5). The digital error signal is sent 
through the servo filter, Fals? which includes a 470 /as 
processing delay. The resultant digital control signal 
is converted back to an analog voltage via a digital-to- 
analog converter (DAC), and the output analog control 
signal is used to actuate on the end test mass via electro- 
magnetic actuators (Atm)- The overall open loop gain 
of this loop is roughly 



Gals — ^dfd ^als ^tm, 



(5) 



since the effect of the AUX-cavity loop and Havx in 
parallel is an overall flat frequency response that does 
not affect the overall open-loop gain. 

When the loop is closed, fluctuations in the frequency 
of the beat note are suppressed by acting on the cav- 
ity length. Any external disturbances that produce fre- 
quency shifts common to both the AUX and PSL beams 
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Fig. 5. Block diagram of the model of control scheme and noise sources for the experiment. The colored blocks 
correspond to the colored blocks in the setup diagram in Figure 1. The various individual control elements are 
described in the text and in appendix A. 



will then be suppressed by the closed loop suppression 
factor, 1/(1 + Gals)- The transfer function between dis- 
turbances common to both wavelengths and residual dis- 
placement can be seen in Figure 8 in appendix A. 

E. PSL PDH sensing 

This block represents the direct out-of-loop measurement 
of the residual displacement noise of the cavity relative 
to the length of the MC. As opposed to the external dis- 
turbances common to the AUX and PSL beams that are 
suppressed by the arm length servo, any external distur- 
bances sensed differentially between the AUX and PSL 
beams will transmit directly to this sensor and contribute 
to any residual displacement noise. The transfer function 
between differential external disturbances and residual 
displacement can be seen in Figure 8 in appendix A. 

There are noise sources here related to readout, such 
as ADC noise, and shot noise and dark noise in the MC 
transmitted photo detection process, but they are found 
to be insignificant relative to other noise sources and are 
therefore omitted. 



4. Performance and Noise Analysis 

The primary usefulness of the multi-color readout is that 
it enables us to precisely adjust the arm length and hold 
it at a desired value independent of the state of the rest 
of the interferometer (as discussed in section 2). Fig- 
ure 6 shows a sweep of the cavity length feedback offset 
through the cavity resonance of the 1064 nm PSL beam. 
The figure demonstrates that the length detuning can 
be cleanly and smoothly brought to zero, at which point 
the 1064 nm beam is fully resonant in the cavity. The 
top plot is the amount of the detuning in terms of the 
beat frequency observed at the fine DFD output. The 
middle plot shows the 1064 nm intracavity power as it 
passes through resonance with the cavity. The bottom 
plot shows the PSL PDH error signal. 

The residual arm displacement measured in the out- 
of-loop PSL PDH error signal has a root mean square 
(RMS) of 23.5 pm, integrated from 1 kHz to 10 mHz. 
The amplitude spectral density and RMS of this resid- 
ual displacement are shown as the solid and dashed red 
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Fig. 6. Sweep of cavity length control offset through the 
1064 nm (PSL) resonance of the cavity. A detuning of 
100 kHz corresponds to a cavity displacement of 6.7 nm. 



curves respectively in Figure 7. The measured RMS falls 
below the aLIGO requirement of 1.3 nm RMS [25], which 
is based on the the line width of the arm cavity for the 
1064 nm wavelength. 

Figure 7 also shows the overall noise budget of the ex- 
periment, i.e. an accounting of all noise sources that are 
thought to affect the performance of the experiment. To 
determine the contribution from a particular source we 
first calculate, estimate, or measure the power spectrum 
of the noise, at its source (designated by a in Fig- 

ure 5), and propagate the amplitude spectrum through 
the control model to produce an amplitude spectrum in 
the "residual displacement" output. The result is 



n(/)=X(/)vW), 



(6) 



where X{f) is the transfer function from the noise source 
to the out-of- loop- measured residual displacement. In 
the rest of this section we will describe the contribution 
from each of the sources shown in Figure 7. 

The total noise spectrum accountable from the budget 
is shown as the solid blue trace in Figure 7. The fact 



that the blue trace lies below the red trace over much of 
the band indicates a discrepancy in the noise accounting. 
The measured displacement spectrum is limited by a 1 //- 
shaped noise at low frequencies, and a white noise above 
100 Hz. It is unknown at this time where these limiting 
factors originate. 

Table 3 in Appendix B describes all of the fundamental 
constants, experimental values, and material properties 
used in this section. 

A. Fundamental cavity noise sources 

This section describes various fundamental noise sources 
in the cavity being measured. While most of these noises 
can't be measured directly, their levels can be estimated 
based on analytical models of the underlying physics. 

1. Seismic noise 

Seismic noise, while dominant across much of the band of 
interest, is suppressed by the cavity length control loop. 
The light gray "unsuppressed" spectrum in Figure 7 is an 
estimate of the free-swinging cavity motion, V^S*^, based 
on the in-loop error signal when the cavity is locked to the 
PSL. This spectrum is expected to be entirely dominated 
by seismic noise at frequencies below 100 Hz, and reaches 
a level of roughly 10~^ m/V^Hz below 1 Hz. The peak 
at 1 Hz is due to the pendulum resonance of the optic's 
suspension system, while the peak at 3 Hz is due to the 
resonance of the vibration isolation stack which supports 
the optical table in the vacuum chamber. 

The dark gray "expected suppressed" trace is the "un- 
suppressed" convolved with the closed-loop transfer func- 
tion from the common cavity displacement input to the 
residual displacement. This represents the expected con- 
tribution of seismic noise to the residual displacement 
once the arm length servo loop is closed. Since this mo- 
tion is common to both laser frequencies its contribu- 
tion is suppressed by the servo to an expected level of 
^ 10~^^ m/v^Hz across most of the band. However, it is 
nonetheless found to be one of the main contributors to 
the resultant arm displacement around 10 Hz. 

2. Suspension thermal noise 

The cylindrical cavity mirrors are suspended from a sin- 
gle wire loop clamped at the top of a suspension cage. 
The length of the pendulum is 25 cm, which leads to a 
fundamental pendulum frequency of fp = 1.0 Hz. Sus- 
pension thermal noise originates from thermal fluctuation 
of these suspension wires. The noise is well modeled [26] 
and its power spectrum is expressed as 
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where ks is the Boltzmann constant, T is the mean tem- 
perature, and / is the frequency. Y is the admittance 
of the suspended mirror due to an external force and is 
described by the transfer function 
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Fig. 7. Residual displacement noise at 1064 nm, and noise budget of the locked cavity. The right vertical axis indicates 
the corresponding frequency noise at 1064 nm. The solid red curve is the overall residual cavity displacement measured 
by the out-of-loop PSL PDH detector, while the red dashed curve is the residual displacement RMS integrated from 
high frequency. The colored dash-dot curves represent the estimated noise contributions from various fundamental 
noise sources, while the solid colored curves are the measured or estimated levels of the various technical noise sources. 



where M is the pendulum mass, fp is the pendulum fre- 
quency, and (j)p is the pendulum loss angle. As with seis- 
mic noise, this noise is common to both AUX and PSL 
beams so its contribution is suppressed by the arm length 
servo. 

The expected suspension thermal noise level in our ex- 
periment is shown as the dot-dashed orange curve in Fig- 
ure 7. The peak at 1 Hz is excess noise due to the pen- 
dulum resonance. The dip at 3.2 Hz is due to a resonant 
gain in the arm length servo used to suppress the contri- 
bution from the primary mode of the vibration isolation 
stack, while the dip at 16.5 Hz is resonant gain used to 
suppress the contribution from the bounce mode of the 
optic suspension. 

3. Coating thermal noise 

Noises associated with thermal fluctuations in the mir- 
ror's high-reflectivity coatings are an important limiting 
noise source in LIGO. While they are not expected to be 
a notable contribution to our result, we touch on them 
here for completeness. 



There are two important coating thermal noise sources: 
Brownian noise comes from thermal vibrations associ- 
ated with mechanical losses in the mirror coating. The 
combined Thermo -refractive and thermo- elastic noises, 
jointly referred to as thermo-optic (TO) noise, affects the 
laser fleld as it interacts with the high-reflective coating. 

Unfortunately, calculating the effect of these noises in 
the presence multiple light wavelengths is not trivial. As 
mentioned in section 3 E, only the noises sensed differen- 
tially between the AUX and PSL beams, and therefore 
not suppressed by the arm length servo, will show up as 
residual displacement noise in our experiment. Calculat- 
ing the differential effect accurately would therefore re- 
quire a fully coherent analysis at both wavelengths, which 
we will not attempt here. Instead we make the very naive 
assumption that the difference between what is sensed by 
the AUX and PSL beams is entirely attributable to the 
difference in their spot sizes on the mirror surfaces. We 
then calculate the differential thermal noise contributions 
based on this differential spot area. 
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The Brownian thermal noise spectrum is given by [27]: radiative transfer to the optic surface, 



ikeT (pes (1 - 
27r/ 



(9) 



where P is the Poisson ratio of the substrate, (/)eff is the 
effective loss angle of the coating, E is the Young's mod- 
ulus of the substrate, and a is the area probed. The 
resultant residual displacement spectrum from coating 
Browning noise in our experiment is shown as the dot- 
dashed brown curve in Figure 7. 

For the thermo-optic (TO) noise contribution, we fol- 
low the coherent treatment proposed in [23]. A Gaussian 
beam illuminating a mirror senses thermal fluctuations 
in the coating resulting in the noise power spectrum: 



(10) 



where k is the thermal conductivity, and s is the heat ca- 
pacity per volume. The overall thermo-optic noise spec- 
trum is then 



5to(/) = S^l{j) Tte (xfsmAad - ~^\y 



(11) 



where Ftc is a correction due to the finite thickness of the 
coating layers, Aa is the difference in effective thermal 
expansion coefficient between the coating and substrate, 
Xfsm is a correction due to the finite mirror size, d is the 
thickness of the layers, ^ is the effective thermo refractive 
coefficient, and A is the beam wavelength. The residual 
displacement from thermo-optic noise is shown as the 
lower dot-dashed blue trace in Figure 7. 

4. Couplings with ambient temperature fluctuations 

Potentially more significant than the inherent thermo op- 
tic noise contribution at low frequencies is the thermo 
optic contribution from low frequency ambient tempera- 
ture fluctuations coupling directly to the mirror coating. 
Thermal fluctuations in the mirror coating due to am- 
bient temperature fluctuations in the lab can be signifi- 
cantly higher than those from thermo optic excitations. 
These fluctuations dominate the thermo optic noise spec- 
trum at low frequencies. 

To estimate the thermo optic noise contribution from 
ambient temperature we start with the same thermo op- 
tic noise description in equation 11. But instead of using 
the thermo optic fluctuations from equation 10, 5'xo ' 
instead use an estimated thermal spectrum given by: 



sf?{.f) = mf) c{f) jiff 



(12) 



where ST{f) is the amplitude spectrum of the ambient 
temperature fluctuations in the lab environment. 



5T{f) = 3 X 10" 



0.01 Hz\ K 



/ 



Hz 



(13) 



C{f) is the transfer function through the vacuum enve- 
lope, described by a single Of.l Hz pole, and j{f) is the 



Jif) 



(14) 



where e is the emissivity of the coating, a is the Stefan- 
Boltzmann constant, T is the mean temperature, k is the 
thermal conductivity of the substrate, p is the density of 
the substrate, and c is the specific heat capacity. The 
contribution from this affect is shown as the pink dot- 
dashed curve in Figure 7. 

B. Technical noise sources 

This section describes the contribution from various tech- 
nical noise sources that can be measured directly in the 
experiment. 

1. Laser frequency noise 

Frequency noise associated with the AUX and PSL lasers 
is generally suppressed by the control loops that keep the 
lasers locked to the main arm and mode cleaner cavities. 
However, since all control loops are coupled together at 
some level, there is a possibility of laser frequency noise 
contributing to the measured residual displacement noise. 

As discussed in section 3B frequency noise from the 
PSL is significantly suppressed and can therefore be ig- 
nored. However, coupling from the AUX laser is at a 
much higher level. We estimate its contribution by ob- 
serving the residual noise in the AUX PDH error signal 
while the AUX-cavity loop is locked. We then assume 
that this noise is due almost entirely to unsuppressed 
laser frequency noise fluctuations. The resultant con- 
tribution from this noise is the dark green "AUX PDH 
error" trace in Figure 7. 

2. Second-harmonic generation noise 

Noise due to the second-harmonic generation process 
is assumed to be added to the frequency noise of the 
frequency-doubled laser beam. The upper limit of the 
noise level is assumed to be 1 x 10~^/ Hz/\/Hz [28]. Laser 
frequency doubling happens in two places in our exper- 
iment: in the AUX laser output and on the PSL beam 
for the beat note detection. The contribution from the 
PSL doubling in the beat note detection is much more 
significant, so it is this level that is shown as the light 
green curve in Figure 7. 

3. Shot noise and detector dark noise 

Both shot noise and dark noise appear as white noise (in 
the detection band) in the broadband RF photodetectors 
used in the experiment. When measuring the frequency 
of a signal, the measured voltage noise can be converted 
to frequency noise on the detected signal by [29]: 



Q f4^\ 26V(/),2 
^PDUj — / : 



(15) 



where 5'v is the input-referred voltage noise, and V^if is 
the voltage amplitude of the main RF signal. 
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The most dominant contribution from these noises 
comes from the beat note detection photodetector. The 
dark current noise level of the photodetector used is 
12 pA/V^ between 10 MHz and 80 MHz. The inci- 
dent power on the PD is 200 fiW which produces 60 /iA 
of DC photocurrent, corresponding to a shot noise level 
of 4 pA/V^Hz. The resultant frequency noise spectra at 
the detector input for these noise sources are shown as 
the purple (dark) and magenta (shot) traces in Figure 7. 

4. Frequency discriminator noise 

The comparator in the delay-line frequency discrimina- 
tor adds white noise during the process of reshaping RF 
signals into square waves. The noise level is measured 
from the output of the DFD while being driven by a 
pure RF sine wave. The level, referred to the input 
of the DFD, was found to have a total contribution of 
]^Q-i4 _ io~^^ m/v^Hz after applying the loop correction 
factor (olive curve in Figure 7). 

5. ADC noise 

ADC noise is easily measured directly by terminating the 
inputs to the analog filters that whiten the signal before 
digitization, and then measuring the spectrum digitally. 
The affect of the of the ADC whitening is compensated 
for within the digital system. The resultant contribution, 
referred to the input of the DFD, sees the same loop 
correction factor as the DFD (cyan trace in Figure 7). 

6. DAC noise 

DAC noise is directly measured by digitally generating a 
3 Hz signal, representing the peak frequency of the error 
signal while locked, and then measuring the output noise 
spectrum. The resultant noise contribution is mostly flat 
at a level of about 1 /iV/A/Hz and is shown as the pale 
blue curve in Figure 7). The dips at 3.2 Hz and 16.5 Hz 
are due to the affect of the resonant gains stages discussed 
in section 4 A 2. 

C. Scaling noise sources for Advanced LIGO 

In this section we look briefly at how various noise sources 
in our experiment can be scaled to Advanced LIGO. We 
find that certain noise contributions will be more promi- 
nent in aLIGO, but that they are addressed in the aLIGO 
design such that they should not pose a significant prob- 
lem. 

1. Frequency noise 

From the relation between frequency noise and length 
fluctuations expressed in equation 2 we can see that 
the 100-times longer arm cavities of aLIGO means that 
aLIGO will be 100 times more sensitive to laser frequency 
fluctuations. For a displacement noise requirement of 
1 nm RMS, the beat note frequency stability require- 
ment goes from 8.8 kHz in the 40m experiment to 83 Hz 
in aLIGO. This puts a much stricter requirement on the 
frequency noise of the AUX laser and SHG noises. 

Advanced LIGO will mitigate this issue in a couple of 
different ways. First, aLIGO will phase lock the AUX 



laser to the PSL frequency through the use of a fiber- 
based phase- locked loop (PLL). This will improve noise 
below 50 Hz, while making it worse at high frequencies. 
The increased high frequency noise can then be addressed 
through optimization of the servo controls. The gain of 
the AUX-cavity loop can be increased to suppress the 
excess noise from the AUX laser, and the bandwidth of 
the arm length stabilization loop can be decreased so that 
any residual noise will not be injected into the cavity 
motion. 

2. Readout and frequency discriminator noise 

Readout and electronics noise sources should become less 
severe in aLIGO. The interferometer response will gen- 
erally grow in proportion to length, resulting in a higher 
SNR against these noise sources. 

The frequency discriminator, on the other hand, gen- 
erally does not scale with the base line length since it 
reads out the frequency of the beat note rather than the 
optical phase. For this reason the readout noise of the 
discriminator will make a 100 times larger contribution 
to the noise budget than it does in our experiment. This 
could likely be the limiting noise source with a frequency 
noise level of 1 Hz/V^Hz at 10 Hz. The situation can be 
improved by using a small-range discriminator such as a 
VCO-based PLL or a longer cable in the DFD. 

3. Seismic length fluctuations 

Length fluctuations due to seismic noise should become 
somewhat easier to handle in aLIGO since the test masses 
will be far more isolated from ground vibration due to 
sophisticated aLIGO seismic isolation systems. Depend- 
ing on how large the residual seismic fluctuations are the 
unity gain frequency of the arm length servo loop should 
be able to be lowered. This is generally good since it 
avoids injection of undesired control noises at high fre- 
quencies. 

5. Future Work 

Besides precise control of suspended Fabry-Perot cavi- 
ties, multi-wavelength readout also has the potential to 
improve performance of other optical systems. Here we 
present futuristic ideas which can potentially reduce the 
fundamental noise sources such as quantum noise and 
mirror thermal noise through the use of multi- wavelength 
readout, as well as an idea to precisely characterize an 
optical cavity. 

A. Multi-wavelength readout for manupulating the 
quantum noise limit 

By resonating multiple laser beams with different wave- 
lengths in a single interferometer, traditional quantum 
noise limits can potentially be modified. 

One example is cancellation of quantum back-action 
for GW detectors [30]. Imagine a main carrier field res- 
onating with high power in an interferometer arm cav- 
ity, and a low-power auxiliary laser beam with a differ- 
ent wavelength resonant only in the interferometer vertex 
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(anti-resonant in the arm cavities). The high-power main 
carrier field would produce quantum radiation pressure 
noise on the test masses. The low-power auxiliary laser 
beam, on the other hand, would sense only the differential 
motion of the two input test masses, and therefore not 
be sensitive to gravitational wave signals. An optimal 
combination of the two carriers beam with Wiener fil- 
ters could then be used to cancel the low-frequency back 
action noise, while not losing information from gravita- 
tional waves. 

Another idea is to resonate both carrier wavelengths in 
the arm cavities. The design of the optics could be made 
such that the optical properties for the different wave- 
lengths are different leading to different frequency sensi- 
tivities for the two beams. For example, the input test 
masses could have higher transmittance for one wave- 
length over the other. By tuning different wavelengths 
and optimally combining their outputs, one may be able 
to shape the quantum noise spectrum in a much more 
flexible way than in the single wavelength case. 

Multiple beams with different wavelengths could also 
potentially be used to manipulate the dynamics of test 
masses in optical cavities. In particular, it is well known 
that the multi-bounce laser fields in optical cavities mod- 
ify the dynamics of the cavity mirrors via radiation pres- 
sure [31]. Under appropriate conditions this can result 
in modifications of the opto-mechanical coupling and a 
higher response against optical phase changes. For ex- 
ample, in gravitational wave detectors that use signal- 
recycling cavities, two wavelengths of light appropriately 
detuned from the resonance of the signal-recycling cav- 
ity can result in a radiation pressure force that reduces 
the effective inertia of the test mass at low frequencies. 
This can significantly amplify the response of the interfer- 
ometer to gravitational- wave signals [32]. It is therefore 
possible that a mult i- wavelength technique could allow 
for surpassing the standard quantum limit over a broad 
frequency band. 

B. Thermal noise estimation 

In a frequency regime where sensitivity is strongly lim- 
ited by mirror thermal noise, it may be possible to use 
two different laser wavelengths to differentially sense the 
thermal noise in the two fields. It may then be possible 
to combine signals from the two lasers in order to yield 
one data stream representing the cavity length fiuctua- 
tions and another with purely the thermal noise. This 
technique may be capable of giving a moderate decrease 
in the effective thermal noise. However, the differential 
frequency noise between the wavelengths would need to 
be investigated more precisely. 

C. Precise cavity mode characterization 

Since multi-wavelength metrology enables us to detune 
the laser frequency from one of the cavity resonances in 
a quasi-static manner, various longitudinal and spacial 
characteristics of the cavity can be precisely inspected. 



Precise scanning of a resonance can provide a measure of 
cavity finesse, while scanning over multiple free spectral 
ranges gives us a measurement of the absolute length of 
the cavity. The frequency spacing of transverse spatial 
modes can be obtained by inspecting the resonances of 
the fundamental and higher-order modes, therefore pro- 
viding information about the cavity geometry and the 
figure error of the cavity optics. The power transmitted 
during scanning can also tell us the mode matching effi- 
ciency between an incident beam and cavity eigenmodes. 

6. Conclusions 

Using multiple lasers, we have demonstrated a tractable 
strategy for sensing the cavity lengths of a complex 
interferometer for gravitational wave detection. This 
method can be replicated and applied to any of the detec- 
tors in the upcoming worldwide network (LIGO, Virgo, 
KAGRA). 

The noise limits are now well understood and well be- 
low the requirements necessary for aLIGO. The mirror 
masses for these new detectors are 30-40 kg and they 
have thermal time constants of several hours. Long pe- 
riods spent without the interferometer locked introduce 
enormous thermal transients in the interferometer which 
perturb the delicate operating state. The use of this tech- 
nique should allow for a significantly higher duty cycle in 
the future. 

In addition, the technique has promise to improve 
the sensitivity of the next generation of interferometers 
through the use of mult i- wavelength readout to partially 
cancel some of the thermal and quantum noises which 
limit the more conventional designs. 

Appendix A: Control Model Transfer Functions 

Table 2 lists all of the zeros, poles, and gains for the 
various control elements in the control model described 
in section 3 and Figure 5. For a system with Pm poles and 
Zn zeros (both specified in Hz) and gain /c, the transfer 
function would be given by: 



X{f) = k 



n„(i+i//zn) 
n^(i+i//p„)' 



(Al) 



Figure 8 shows Bode plots (amplitude and phase) of a 
couple of the key transfer functions from the full control 
model. 

In addition to the blocks represented in Table 2, the 
blocks labeled "L i^" in Figure 5 represent the conver- 
sion from displacement {dL) to frequency {dv) described 
in equation 2, i.e.: 



(A2) 



where v = c/A. 
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element 


zeros (Hz) 


poles (Hz) 


gain 


Caux 


- 


18.5k 


1 


^AUX 


- 


- 


5.0e-6 


i^AUX 


1.0, 100, 10k 


0.1m, 1.2, 2.0 


2.1e8 


^AUX 


- 


100k 


5.0e6 


^DFD 


- 


- 


2.16e-7 




40, 40, 


1.0m, 500, 




i^ALS 


1.655 ± 2.739i. 


0.052 ± 3.200i, 


1.0e6 




3.320 ± 16. 163i 


0.105 =b 16.500i 




^TM 


- 


0.1±0.995i 


8.0e7 


C'pSL 


- 


3.8k 


1 


^PSL 






2.5e-5 


FpSL 


4k, 20k, 20k 


40, Ik, Ik 


2.3e4 


^PSL 




100k 


5.0e6 


SHG 






2 



Table 2. Zeros, poles and gain of the control model 
blocks. 



ALS open loop, measured 
-ALS open loop, model 
-common residual 
-differential residual 




symbol 


name 


value 


SI unit 


ks 


Boltzmann's constant 


1.38e-23 




a 


Stefan-Boltzmann constant 


5.67e-8 


W m-^K-^ 


T 


mean temperature 


290 


K 


A 


beam wavelength 


1.064e-6 


m 


M 


cavity mirror pendulum 
mass 


0.243 


kg 


fp 


pendulum frequency 


1.0 


Hz 




pendulum loss angle 


1.7e-4 


_ 


p 


Poisson ratio of (substrate) 


0.167 




E 


Young's modulus (sub- 
strate) 


7.27el0 


N m-^ 


K 


thermal conductivity 


1.38 


W m-^ K-^ 


S 


heat capacity per volume 


1.62e6 


m-3 


c 


specific heat capacity 


740 


J kg-^ K-^ 


e 


emissivity 


0.9 




P 


density 


2202 


kg m~^ 


ITM 




effective coating loss angle 


8.65e-8 




Aa 


effective thermal expansion 
difference 


3.59e-6 






effective thermal refraction 


2.35e-6 




ETM 


(/>eff 


effective coating loss angle 


1.24e-7 




Aa 


effective thermal expansion 
difference 


4.54e-6 






effective thermal refraction 


l.lle-6 





10 10 
Frequency [Hz] 



Table 3. Values of fundamental constants, and material 
properties and variables for the 40m mirror coatings. If 
not specified, material properties are for the mirror coat- 
ing. 
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Fig. 8. Bode plot of control model transfer functions. 
The green trace is the full open-loop transfer function of References 
the arm length servo control loop. The blue and red 
traces are the transfer functions from external distur- 
bances to residual displacement, for AUX/PSL common 
and differential sensing. 
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Appendix B: Symbol Definitions and Values 2. 

Table 3 shows the values used for various variables in the 
text, including all fundamental constants, experimental 
values, and material properties. 
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